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Abstract 

We find exact power-law solutions for scalar-tensor theories and clarify the conditions under which they 
can account for an accelerated expansion of the Universe. These solutions have the property that the signs 
of both the Hubble rate and the deceleration parameter in the Jordan frame may be different from the signs 
of their Einstein-frame counterparts. For special parameter combinations we identify these solutions with 
asymptotic attractors that have been obtained in the literature through dynamical-system analysis. We 
establish an effective general-relativistic description for which the geometrical equivalent of dark energy is 
associated with a time dependent equation of state. The present value of the latter is consistent with the 
observed cosmological "constant" . We demonstrate that this type of power-law solutions for accelerated 
expansion cannot be realized in f(R) theories. 



* E-mail: cedumagalhaes@hotmail.com 
t E-mail: winfried.zimdahl@pq.cnpq.br 



1 



I. INTRODUCTION 



The past decade has seen a continuously growing activity, both theoretical and observational, to 
achieve a consistent picture of the large-scale cosmological dynamics. A huge amount of data has 
been accumulated which directly or indirectly seem to back up the conclusion, first obtained in 
that our current Universe entered a phase of accelerated expansion. Direct support is provided by 
the luminosity-distance data of supernovae of type la [2] (but see also [3i]), indirect support comes 
from the anisotropy spectrum of the cosmic microwave background radiation (CMBR) [4], from 
large-scale-structure data [5], from the integrated Sachs- Wolfe effect [6(, from baryonic acoustic 
oscillations [?J and from gravitational lensing 

A satisfactory explanation of this originally surprising result has not yet been achieved so far. 
The component responsible for the accelerated expansion is called dark energy (DE) but its physical 
nature remains unclear. The first guess, an effective cosmological constant, remained the favored 
option until today and gave rise to the ACDM model which also plays the role of a standard 
reference for any research in the field. According to the interpretation of the data within this 
model, our Universe is dynamically dominated by a cosmological constant A which contributes 
roughly 75% to the total cosmic energy budget. Roughly 20% are contributed by cold dark matter 
(CDM) and only about 5% are in the form of conventional, baryonic matter. Because of the 
cosmological constant problem in its different facets, including the coincidence problem (see, e.g., 
0, [l(3] ) , a great deal of work was devoted to alternative approaches in which a similar dynamics as 
that of the ACDM model is reproduced with a time varying cosmological term, i.e., the cosmological 
constant is dynamized. An important class here are the quintessence models which, in a sense, 
mimic scalar field inflationary models for the early universe, albeit on a very different energy scale. 
Also scalar fields with a non-standard kinetic term, called k-essence, have attracted attention. A 
further class are unified models of the dark sector, i.e. models, for which a single component plays 
the role of CDM in the past and the role of a mixture of CDM and DE at the present time. The 
generalized Chaplygin gas is t 
models belong to this category 



re prototype for this kind of models [111, [12j. Also viscous fluid 
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14] . Common to all these approaches is that Einstein's General 



Relativity (GR) is assumed to be the valid gravitational theory up to the largest cosmological 
scales. But the apparently strange nature of DE and the fact that both dark matter (DM) and DE 
manifest themselves only gravitationally has also provoked the exploration of alternative theories 
of gravity. A major line of investigation in this context are scalar-tensor theories. Brans-Dicke 
theory, based on ideas of Mach and Jordan [l5 ] , is the prototype of a scalar tensor theory [u], [l^] • 
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According to Mach's principle, the gravitational coupling depends on the mass distribution in 
the Universe. Hence the (effective) gravitational "constant" may be a variable quantity since 
the mass distribution may vary. Based on this motivation, Brans and Dicke introduced a non- 
minimally coupled scalar field in order to modify the Newtonian gravitational interaction. Also 
more general scalar-tensor theories that have been developed subsequently, are characterized by 
the fact that the gravitational interaction is mediated both by a metric tensor and a scalar field. 
The interest in modified theories of gravity is triggered by the expectation to describe the late- 
time accelerated expansion of the Universe without a DE component 1814201] . Instead, it is the 



geometrical sector which is supposed to provide the desired dynamics which is then gravitationally 
induced 
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- 241] . This class of approaches can be seen as a geometrization of DE. Fundamental 
quantum theories involving extra dimension typically seem to predict 4-dimensional scalar-tensor 
theories 



251 ] . Modified gravitational theories in which the Ricci scalar R in the Einstein-Hilbert 
action for GR is replaced by a function of R are called f(R) theories. This line of investigation 
has attracted particular attention as an alternative to GR ([261429]). f(R) theories were shown to 
be special cases of scalar-tensor theories With a suitable self-interacting potential and 

an appropriate coupling to matter fields (chameleon effect) these theories are potential candidates 



for a geometrical description of DE 



33H37I] (see also 



38, 



viability of f(R) dark energy models were formulated in 



these theories see 
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39l|). Conditions for the cosmological 



40]. For recent reviews on the status of 



Scalar-tensor theories are formulated both in the Einstein frame and in the Jordan frame, the 
two being related by conformal transformations. They are considered to be physically equivalent, 
although there may occur differences concerning the General Relativity (GR) limit [44]. On the 
quantum level the equivalence may be lost. For a discussion of apparent interpretation problems 
see 4514471] . Various asp ects of scalar-tensor theories in general or subclasses of them have been 



investigated 
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4814521]. Scalar-tensor theories have also been used in attempts to provide a 



geometrical explanation of DM 



23, 



29, 



5j455(]. 



Scalar-tensor theories are more complex than GR. Even if the symmetries of the cosmological 
principle are imposed, scalar-tensor theories do not admit simple solutions that could be compared 
with, say, the GR power-law solutions for the cosmic scale factor for fluids with constant equation 
of state (EoS) parameters. Investigations in the literature frequently rely on a dynamical system 



5J460I]. 



analysis with the aim to find critical points, equivalent to asymptotic power-law solutions 
But to the best of our knowledge, exact power-law solution do only exist for special classes of 
"curvature quintessence" and without a matter component 6l[ but not in the general case. It is 
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the aim of this paper to partially fill this gap and to derive simple exact solutions which allow 
for a transparent discussion of at least some of the aspects of these theories. To this purpose 
we start looking for solutions with a constant ratio of the energy densities of the matter and the 
scalar field components in the Einstein frame. The resulting scaling solutions imply a relation 
between the equation of state (EoS) parameters of the components and the interaction-strengt 



parameter. For previous studies of scaling solutions of the cosmological dynamics see, e.g., 
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68]. The Einstein- frame solutions are then transformed into power-law solutions of the Jordan 
frame, representing a one-component description of the dynamics that can be confronted with 
results from dynamical system analysis. In a next step we interpret the Jordan- frame solutions as 
solutions of an effective Friedmann equation within GR. Separating a conserved matter part from 
the total substratum enables us to identify the effective EoS of the (conserved) remaining part that 
now is regarded as the equivalent of DE. This two-component description is then compared with 
the ACDM model. 

The paper is organized as follows. Section HT1 recalls the basic dynamics of scalar-tensor theories 
both in the Jordan and in the Einstein frames. In Section IIIII we obtain scaling solutions in the 
Einstein frame. The corresponding Jordan-frame solutions are found in Section IIVI Based on 
these results we discuss the effective EoS of the cosmic medium and point out the relation of our 
approach to the ACDM model in Section [V] In Section IVT1 we test the obtained power-law solution 
against supernova-type-la data. Section IVIII provides a summary of the paper. 

II. BASIC DYNAMICS 

We start by reviewing the basic relations for scalar-tensor theories where we adopt the notation 

of 



58J. Scalar-tensor theories are based on the (Jordan- frame) Brans-Dicke type action 
S = J d A x^—g 

+ / d A Xy/^gL m {g^ u ) , (1) 



(*) R - ^~^)_}_ (d_F\ 2 2 _ y 

2k 2 K ' 16k 2 /3 2 F($)\d$J { } { } 



where k 2 = 8irG and (3 is a coupling constant. The quantities g and R are the determinant and the 
curvature scalar of the metric tensor g^ u , respectively. L m denotes the matter Lagrangian. F > 
is a function of the scalar field with a potential V(3>). With the help of the transformations 



9tiv 



e 



-a § ; = g 2 ^^ , = F($) 2 V{<j>) , (2) 
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one obtains the Einstein frame action 
S = I d^x\J ' —g 



(3) 



An explicit expression for is required to specify $ as a function of 4>. Throughout the paper, 

quantities with a tilde refer to the Einstein frame, quantities without tilde have their meaning in 
the Jordan frame. 

In the following we restrict our attention to the application of the general theory to homogeneous 
and isotropic cosmological models with flat spatial sections. Moreover, we assume that the matter 
part can be modeled by a perfect fluid. Under these conditions the relevant equations are 



H 2 



8nG , . HdF 1-4/3 2 1 / dF N 



3F 



[Pm+vm- 



+ 



- 



F dt 16/3 2 F 2 \ dt 



where H = is the Jordan- frame Hubble rate, 



a dt 

dH 



1 _ „, 1 d 2 F 1H dF 1-4/3 2 1 / dF\ 2 



+ 



— 



dt ~ 2F 87TG{Pm+Pm) 2F dt 2 ' 2Fdt " 16(3 2 F 2 \dt) ' 



-3- 



(4) 



(5) 



and 



AV{<S>) - 2V> + (p m - 3p m ) 



dp r . 



dt 



+ 3H ( Pm +p m ) = 



The corresponding Einstein-frame relations are 



H z = — 



with the Einstein- frame Hubble rate H = i ^ , 



dH k 2 



+ Pm 



dt 



(6) 



(7) 



(8) 



(9) 



and 



d 2 6 ~ dd> ~ 



+ 3H (p m + p n 

dt 



K/3 (p m - 3p m ) , 



3 K/3 ~dl ^ m ~ 3 ^ 



(10) 



(11) 
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The time coordinates and the scale factors of the Robertson- Walker metrics of both frames are 
related by 

dt = e-V^^ d t and a = e~ VW»*a , (12) 
respectively. The matter quantities transform into each other via 

Pm — e Pm j — t- fm ■ V 10 7 

This means 2s. = 2™. i. e the EoS parameter remains invariant. 

III. A TWO-COMPONENT DESCRIPTION 
A. General relations 

The set of Einstein-frame equations ([8]) - (jlip suggests an effective two-component structure in 
which matter interacts with a scalar field. One may attribute an effective energy density p^ and 
an effective pressure to the scalar field by 

1 fdd)\ 2 ~ , „ ! /.7 n 2 



+ y md w=2 li) - v ' (14) 



respectively. Equations (|TT]) and (|T0|) can be written as 



' ' // '' n + 3tf (1 + «) m ) = -t/f (1 - 3* m ) (15) 



p m dt \ 3 dt 

and 



_L^| + 3H (1 + fy) = (1 - 3^ m ) r , (16) 

p^ dt V o £tt 

respectively, where we have introduced the matter EoS parameter w m = and the Einstein-frame 
EoS parameter for the scalar field = lp. Furthermore, we have defined the ratio of the energy 
densities 

Pm f-in\ 

r = -=~. (17) 

P<f> 

The total energy density and the total pressure are 

P = Pm + P^,, P = Pm+P<f> => ^ + 3H(p + p)=0. (18) 



The set of equations (|15p and (|16p is reminiscent of interacting quintessence models (see, e.g., 



681), 



in which the scalar field is part of the energy-momentum tensor within standard GR and interacts 



with the matter component in a specific way. From this point of view (|15[) and (|16p can be 
regarded as a pure GR-based model. Then the parameter r describes the ratio of dark matter to 
quintessential dark energy. However, such type of interpretation masks the circumstance that in 
the present context the scalar field is a gravitational degree of freedom. But the formal equivalence 
between ([To]) and (fTUj) and corresponding equations in GR may be used to apply solution techniques 
of the latter to the former situation. This will exactly be our strategy in the following subsection. 
A physical interpretation will be given subsequently within the Jordan frame. 

We mention that the coupling parameter /3 is related to the Brans-Dicke parameter ujbd by 



3 1 - 4/3 2 
wbd = - 



1 



(19) 



/3 2 4 2co BD + 3 

It is known that f(R) theories can be regarded as a subclass of scalar-tensor theories, corresponds to 



the special case j3 2 = | in Eq . ([Tl) (cf . 



ruled out observationally 



30j-|32j), equivalent to QJrd = 0. These theories seemed to be 



3-37 



42 



43J , a suitable effective potential may 



3d . [52| , but according to 

result in a sufficiently heavy mass of the scalar field in regions of high matter density (chameleon 
mechanism), so that conflicts with solar system constraints can be avoided. The circumstance that 
scalar-tensor theories are related to a non-linear Lagrangian by a conformal transformation is also 
known as Bicknell-theorem 69], Q] . 



B. Scaling solutions 



From now on we shall focus on the subclass of solutions for the system (|15p and (|16p that admit 
a constant ratio of the energy densities of both components. Corresponding solutions in GR are of 
interest in connection with the coincidence problem, i.e. the question, why the densities of DM and 
DE are of the same order just at the present time. In the present context such connection is less 
obvious since for an adequate interpretation one has also to consider the Jordan-frame solutions 
which will be found below in section HVl The general dynamics of the energy density ratio r is 
obtained by differentiating (fTT|) and using the balances (fT5|) and ([To]) : 



dr 

-z = 3Hr 
dt 



w<f> - w m 



2 „ 1 d6 . 



(20) 



Scaling solutions, i.e. solutions with a constant value of the Einstein-frame energy density ratio r, 
are then characterized by 

dr , 12 



dt 







3 «/3^(l-3tD m ) = 3ff 



W<$> - Wm 

l + r 



(21) 
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Notice that w m = | corresponds to the interaction- free case in (fT5|) and (fT6j) . Under this condition 
a scaling solution can only exist if w$ = | as well. The left-hand side of the second equation (|2"Tj) 
determines the source (loss) terms on the right-hand sides of the balance equations (fT5|) and (fT6|) . 
Consequently, the latter equations become 



;/ . + 3H (ft + ft) = 3ff rft . (22) 

and 

% + 3tf (ft, + vm) = -3H ^LZ^IH p m , (23) 
dt 1 + r 

respectively. Hence, for scaling solutions the coupling is completely specified. No free parameter 
occurs on the right-hand sides of Eqs. (j22[) and (|23f) . The direction of the energy flow depends on 
the sign of zft — w m . For — w m > we have a flux from the m component to the component. 
For vjfj, — w m < it is the opposite (always assuming H > 0). 

While the set of equations (|15p and (|16|) is well known in the literature, the configuration (|22|) 
and (|23|) does not seem to have attracted attention so far. Studying the system (f22|) and (f23j) . 
exploring its consequences for the Jordan frame and discussing implications for the corresponding 
cosmological dynamics are the main aims of this paper. 

For constant values of w m and ?ft equations (f22l) and (l23j) have the solutions 

ft oc p m oc a 1+ 7 . (24) 



Here we have used that integrating Eq. (|2ip yields ev 2 / 3K ^ oc a (i-3<z> m )(i+>-) _ The Friedmann 
equation ([8]) takes the form 3H 2 = 8nGp = 8nG (1 + r) ft which results in 

5a |ll+VFfej (25) 

for the scale factor. Obviously, the correct limits for r — > and r — > oo are consistently recovered. 
The dynamics is that of a substratum with an effective EoS 

= ™0 + ™>m ^ = _2 1 

1 + r ' 3(1 + w cS ) t V 7 

There exists also a contracting solution 

. _ ft f -t\ 3(i+kfr) ~ 2 1 

a = ai J- — =r- , H = , B , (27) 

where ti,t < tf. With (I26p (or (I27p ) the Einstein- frame dynamics is completely solved for the given 
configuration. The similarity to GR power-law solutions is obvious. But as already mentioned, we 



consider the solutions (|26|) here as an intermediate result and postpone a physical discussion to 
the next sections. 

For the deceleration parameter q = — ^p-^f = ~^ ~ Jfi^jr we obtain q = ~ (1 + 3w e g). The 
condition for accelerated expansion in the Einstein frame is 



1 1 , 

w cS < -- 44> w^ + rwm < -- (1 +r) . 



(28) 



3 * 3 

An accelerated expansion can be obtained for w m = and < — | (1 + r). The special case of 



exponential expansion is characterized by 



WcS 



T <^ 1 + + r (1 + w m ) = 44> Wff, + rw m = — 1 — r . 



(29) 



Notice that it^ > —1. Consequently, 



-1 



1 + Uty 



(30) 



i.e., w m < — 1, the component m has to be of the phantom type. A standard matter type solution 
p oc a~ 3 is obtained for u) e g = 0. According to (|26|) . this is realized for any combination w^+ruim = 
0. 

From Eqs. ()21[) and Friedmann's equation it follows that 



1 
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(w m - w<f,y 



2 \ dt J 4 /J2 (1 _ 3 u, m ) z (1 + r) 



(31) 



The potential is given by 
V(6) -. 



1 

2 V dt 



1-?. 



P<p 



(32) 



4/32(l-3u; m ) 2 (l + r)_ 

Consistency between (|32p and V" = i (1 — tt)^) implies the following expression for the interaction 



parameter, 



I3 [ 



2 _ 9 (w m - w^) 2 1 



1 



(33) 



2 (1 - 3w m ) 2 1 + r I + W4' 
As already mentioned, the interaction constant in the balances (|15p and (|16|) is not a free parameter 
for the given configuration, but determined by the equation of state parameters and the ratio r. 
Eliminating f3 via (|33p . the potential may be written as 



V{4>) = Vi exp{-v / 247G 1 



1 + r 



(l + w eS ) {4>-4>i)} 



(34) 
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With (|M|) and (j21j) the scalar-field equation (|T0j) is identically satisfied. Finally we notice that 
combination of (fT9"j) and (f3"3"|) results in 

3 1 (1 - 3w m f „ 

wbd = -- + — — — 2 (1 + r) (1 + u^) (35) 

2 12 K - uty) 

'or the Brans-Dicke parameter. Accelerating solutions for dust with uibd < have been studied in 

This concludes our consideration of scaling solutions in the Einstein frame. In the following 
subsection we shall use the transformations (|l'2p to obtain the corresponding Jordan-frame dynam- 
ics. While this dynamics will also be characterized by power-law solutions, it is not associated with 
a constant ratio of the energy densities of the dynamically relevant components. 



IV. JORDAN FRAME SOLUTIONS 
A. General power-law structure 

The transformations from Einstein's to Jordan's frame are mediated by exponentials of 4>. Via 
relations ()12p one finds 

dt dt V 1 H dt y J 



The Hubble rates are related by 



1 - 3w cS _ -3 - 

H oc a (i-3»m) H . (37) 

1 - 3w m 



This means, the Hubble rates of both frames do not necessarily have the same sign. For w m = 
and w e g > ^, e.g., an expanding solution in the Jordan frame corresponds to a contracting solution 
in the Einstein frame. Similar properties of conformally related frames have been discussed in 
string-theory based pre-big bang scenarios [3]. The explicit relations between the scale factors a 
and a and the time coordinates t and t are 

l-3w cS l-a eff -a m (l+3a cff ) 

a oca 1 - 3 *™ and t oc t a-^ m ){i+a eS ) ^ ( 38 ) 

respectively. These relations encode the correspondence between Einstein frame dynamics and 
Jordan frame dynamics for our power-law solutions. Combination with ()26[) provides us with the 
result 

a oc t ^^es) W ith w e g = . (39) 

1 - 3w cS 
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With (|39p the Jordan-frame dynamics is solved as well. By direct calculation one checks that 
with the Jordan-frame solution (|39p the set of equations ([!])- (HD is satisfied. The solution is 
of the power-law type, corresponding to a constant effective EoS parameter. It is w e g which 
physically characterizes the EoS of the cosmic substratum. Obviously, the parameter w e g differs 
from the Einstein- frame parameter w e g, the letter being an auxiliary quantity. While constant EoS 
parameters certainly cannot account for a continuous transition between different epochs of the 
cosmological evolution, such as the transition from matter to dark-energy type dominance, they 
are useful as simple exact solutions which are valid piecewise. One may associate a conserved total 
effective energy density 

p eS oc a - 3 ( 1+ ^«) (40) 

to the solution ([39]) . The introduction of this quantity allows us to understand our solution alter- 
natively as a GR solution of an effective Friedmann equation 3H 2 = 8iTGp c g for a medium with 
EoS parameter w e f[. We shall come back this point in section PVT. 
The Jordan-frame deceleration parameter can be written as 

" = - 9 = I (1 + 3w M ) , (41) 



(S 2 2 

The relation to the Einstein-frame deceleration parameter is 

_ (1 - 3w m ) (1 + r) _ 1 — 3w m 

q = q-, — ~ — : — t, — ^r-r = q-, — • ( 42 ) 

1 - 3w^ + r (1 - 3w m ) 1 - 3w eS 

Notice that also q and q do not necessarily have the same sign. For the already mentioned case 
w m = and w e ff > | deceleration in the Einstein frame corresponds to acceleration in the Jordan 
frame. For Einstein-frame exponential expansion (|29|) it follows that 

q = -- (1 - 3w m ) , (w^ + rw m = -1 - r) . (43) 

For the degenerate case w m = — 1 we obtain q = — 1 as expected. 



B. Solutions for non-relativistic matter w m = 

For the particularly interesting case w m = (recall that w m = w m ), the second relation (|39l) 
reduces to 

2 ^eff , AA s 

Weft = r. ^-r— ■ (44) 

1 - 3w e g 
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Given the EoS of pressureless matter for the material content of the Universe, Eq. (|44p defines the 
potentially possible physical EoS parameters for the total cosmological dynamics in terms of the 
Einstein-frame parameters. The relation between different intervals of w c g and w e s is shown in 
Tab. HI Interesting special cases in the Jordan frame are: 

(i) a matter dominated phase w e g = 0, realized for w e g = 0. Obviously, this solution requires either 
Wfj, = as well or r — > oo (cf. relation (|26p ). This phase corresponds to the limit between the 
ranges III and IV in Tab. HI According to (|33p . this matter dominated solution has f3 2 = 0. 

(ii) a phase of exponential expansion, realized for w e s = — 1 in (|44p . This requires 

w cS = 1 =4> Wcj, = 1 + r (45) 

and corresponds to the limit between I and II in Tab. HI We emphasize that uyj = — 1 implies 
w e g = +1. While the expansion in the Jordan frame is exponentially accelerated, we have a 
decelerated contraction in the Einstein frame. Recall (cf. Eq. (|37p ) that for w e {[ > | the Hubble 
rates in both frames have different signs. We notice also, that uu > 1 implies a negative potential 

(hi) the onset of accelerated expansion w e g = — | . For this case at the limit between the intervals 
IV and V in Tab. HI we have w e g = — | as well, with = — | (1 + r). 

(iv) an intermediate case, belonging to range I in Tab. HI is w e g = — | , corresponding to w e g = 3 
with ibtf) = 3 (1 + r). 

(v) a phantom EoS w e g = — |, part of range II in Tab. [H This requires w c g = |. 

The existence of power-law solutions with w e g < — | in a universe hlled with pressureless matter 
is possibly the simplest demonstration for the capability of scalar-tensor theories to account for an 
accelerated expansion without a dark-energy component. To the best of our knowledge a power-law 
behavior is known so far only asymptotically but not as an exact solution of the full theory. 

In a next step we investigate, which EoS parameters w e s in (|44p are admitted for (3 2 = |, i.e., 
for /(i?)-type theories. With w m = we find from (J33|) and (f26l) that /3 2 = \ is realized for 



W eff = o 1 ■ (46) 



(vi) /3 2 = j and r = 0. The solutions are w e ff w 2.5 and w e ff « —0.26. Accelerated expansion is 
impossible under these conditions. 

(vii) /3 2 = | and r — > oo. For very large values of r (r > 10 3 ) the solutions of (|46|) are w e // ~ ^ 
and w e ff « 0. The latter reproduces a matter era with a oc 
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591]. For 



(viii) /3 2 = \ and r = 8. The upper sign of f4*6|) yields w e ff = |, corresponding to a oc t 1 / 2 . This 
solution has played a role in the discussion about the cosmic viability of f(R) theories 
the lower sign it follows that w e f f = — ^ . Consequently, there do not exist scaling solutions of the 
type discussed here in f(R) theories that can describe an accelerated expansion of the Universe. 



C. Solutions for w m 7^ 

So far we have considered special cases with w m = 0. But it is also possible to realize equations 
of state of interest under the conditions w m ^ and = 0. For = the general EoS parameter 
in (j39|) specifies to 

r(l-3w m )-l , . 

™eff = Wm—r, o TTT • ( 47 ) 

r (1 - 3w m ) + 1 

This leads to the following quadratic equation for w m , 

9 1 1 — r (1 + 3w P f f) 1 + r . 

wl + o " — W m + -Z-W eff = . 48 

3 r 3r 

We have 

(ix) a matter- dominated universe uyr = 0. There are two solution, w m = and w m = 3^7"- F° r 
r>l the latter approaches w m = 3, i.e., the EoS for radiation. 

(x) exponentially accelerated expansion w e a = — 1. In this case the EoS parameter of the matter 
component is given by 



1 l + 2r 

Wm = -77 

o r 



/ 1 + r 

l± 4 /l + 12r- 



(49) 



(l + 2r)' 

For r> 1 the solutions are w m = —1 and w m = |. For r = 1 one has w m = —1.46 and w m = 0.45. 
For the solutions with w m < it can be argued that now the <j) component describes the (non- 
relativistic) matter and the m component plays the role of dark energy. Different to the cases in 
which the dark energy is geometrized, now it is the dark matter. But again, there exist solutions for 
an accelerated expansion which do not require negative EoS parameters of the matter component, 
i.e., solutions with w c s = — 1 that have = and w m > 0. 

D. Relation to other work 

In a next step we clarify the relation of our solution (|39|) with results from the analysis of 
dynamical systems. The latter typically determines critical points of the cosmic evolution, which 
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EoS Einstein frame 


EoS Jordan frame 


I 


OO > W c ff > 1 


-f > w cS > -1 


II 


1 > w off > | 


— 1 > W c f{ > — OO 


III 


§ > *cff > 


00 > w c ff > 


IV 


> w eS > -i 


> W c ff > -§ 


V 


-3 > w eS > -1 


-3 > W efI > 



TABLE I: Correspondence between different ranges for the Einstein-frame EoS parameter w c ff and the 
corresponding EoS parameter w e e in the Jordan frame for w m = 0. 



5614601] ) . It is therefore 



(in some studies) amounts to finding asymptotic power law solutions (see 
interesting to compare the asymptotic power-law behavior obtained within the context of dynamical 
systems with our exact power-law solutions. As an example we show the consistency of our solution 
with several of the critical points found in 58]. 

0. It 



Point PI in [58( has w e s = 
corresponds to our case (ii) above with r = 0. 



1 and a vanishing fractional matter contribution, — 



Pm 
Pcii 



9-4/3 2 



Point P2 in [5^ has u> eff = and p ^ {3 _^ 2) 

dominated phase of case (i) above. But for (3 2 = | we can also reproduce w e g 
w m = r = and w e s = = 1, corresponding to case (ii) above. 



For /3 2 = it describes the matter 
-1, realized for 



Point P3 in 



58( has w e ff 



3+6/3 and pZ 



0. For our solution 



for f3 with w Tl 



and 



w e & = Wrf) this can be written 



3^3^/2 



Wcff 



3±9\/2- 



(50) 

^ "I, 



For w c s = 1 there is a solution w c s = (upper signs) and another one (lower signs), w e s 
which again reduces to our previous case (ii). 

Finally, Point P7 in [58] has w e s = frff an d = 0. Here, w c s = — 1 can be obtained for j3 = |. 
Again, this requires = 1, i.e., a vanishing potential term. We conclude that our solution (I39p . 
in particular the special case (|44p . is consistent with results from dynamical system analysis. 



V. THE COSMIC MEDIUM 



Now we come back to the effective Friedmann equation with an energy density (|40p . Obviously, 
is a solution of Friedmann's equation 3H 2 = 8irGp c s- The energy density (|4*(J|) . on the other 
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hand, can be seen as the solution of a conservation equation 

dp c s 



dt 



with an effective pressure 



+ 3H(p eS +p eS ) = (51) 



PeS = WeSPeS • (52) 



This set of equations represents a GR equivalent for the power-law solution of the scalar-tensor 
theory. The cosmic medium as a whole is described by an EoS parameter, given by the second 
relation in (]39p . Moreover, we know that a separately conserved matter component with an EoS 
parameter w m is part of the cosmic medium. By direct calculation one confirms that upon using (| 1 3|) 
together with the solution (|24|) the conservation relation ([7]) is consistently satisfied, corresponding 
to 

p m oc a- 3 ^ 1+w ^ . (53) 

Let us define an energy density p x as the difference 

Px = p e a - Pm • (54) 

This means, we assume the total effective energy density p e g to be the sum of a matter contribution 
and this newly introduced and so far unknown x-component. Now we know explicitly both the 
total cosmological dynamics, given by w e s and the dynamics of the separately conserved matter 
subsystem, given by w m . This knowledge allows us to determine the dynamics of the subsystem 
with the energy density p x . Assuming for this component an equation of state p x = w x p x , it is 
then possible to calculate the parameter w x . With p e g = p m + p x and p e g = p m + p x we have 

w eS = w m — +w x — =5- w x = (1 + r eff ) w eS - w m r eS , (55) 

PeS PeS 

where r e g = ^p- is the ratio of the energy densities of both components. Using (|54f) . we obtain for 
this quantity 

pm 

r eS = P™ = . (56) 

Different from the earlier introduced constant energy density ratio r in the Einstein frame, the 
generally time dependent r g is considered to be the ratio of the energy densities of the dynamically 
relevant components within the effective GR description. With the solutions ([53]) and (f4T)|) the ratio 
scales as 

Pcff' 

Prn = PmO Q -3(w m -w rff ) ^ 
Peff PeiTO 
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where is the present ratio of matter energy to total energy. Consequently, for w m = 0, 

^eff _ 1 2q-l 1 2q-l 

WX 1 " 7%To a3Wca 3 1 " ^a^ 1 ' 3 1 _ n m0 (1 + zf~^ ' 1 j 

where z is the redshift parameter and f2 m o = 

As an immediate consequence we find that in the matter era w e g = one has w x = as well. 



7|]. But w x is 



This behavior is reminiscent of specific interacting holographic dark energy models 
not constant for w e s ^ 0. The EoS parameter today is 

Wxo = -r—^. — • (59) 

For the present accelerated expansion, i.e., w e g < one has \w x o\ > \uu e g\, i.e., w x q is more 
negative than w e g, since £l m o < 1- This means, the x-component, here of geometrical origin, 
effectively behaves as dark energy. Now, it is straightforward to make contact with the ACDM 
model. In the latter we have for the total EoS parameter wlcdM) 

wlcdm = 7^ • (60) 

PA + PmO 

With pa = — pa and, observationally, p m o w \pk : the total effective EoS parameter of the ACDM 
model at the present time is wlcdm ~ — |- Identifying the latter tentatively with w e s in ([35]) . we 
find w x o ~ — 1. Consequently, the scalar-tensor-theory solution for the present epoch is compatible 
with the general-relativity based ACDM model, but the equivalent of the EoS parameter for the 
dark energy is time varying. For a future attractor solution w e g = —1, one has w x < —1. For 
a S> 1 the parameter w x approaches w e g from the phantom side. While we have a power-law 
solution for the cosmic medium as a whole, corresponding to a constant EoS parameter w e f[, the 
EoS parameter of the geometrized dark energy is time dependent according to (J35J). The time 
dependence itself is again governed by w e s. 

VI. OBSERVATIONS 

As a final step we test our power-law solution against the supernovae type la (SNIa) data 
and compare the results with the corresponding analysis within the ACDM model. The relevant 
quantity here is the luminosity distance which, in a spatially flat universe, is given by 

d L = (l + z)f^. (61) 

According to (139)) the Hubble rate is 

H = H a~K 1+w ^ =H (1 + z)^ 1+Wca) . (62) 
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Then, taking into account (|4ip . the result of the integration is 

1 + z 1 



d L 



1 



(63) 



H (-g ) 

i.e., the luminosity distance is exactly known. We are interested in applying this expression around 
the present epoch. Expanding in powers of z, we obtain for the luminosity distance up to the third 
order in z 



Hn 



l + \{l-qo)z+ l -{ql-l)z 2 + .. 



(64) 



The corresponding expression of the ACDM model is (cf. Q H) 



>ACDM _ 

L ~ H 



1 + I (1 - q ) z + 1 (3 (q + l) 2 - 5 (q + 1)) z 2 + 



(65) 



Up to second order, the luminosity distances of both models coincide. But they differ at third 
order. 

In the following we will use the sample of 182 SNe la of the Gold06 data set [r\ . The crucial 
quantity for our analysis is the moduli distance fj,, which is obtained from the luminosity distance 
dL by 



fj, = 5 log 



d L 
Mpc 



+25 



(66) 



In order to compare the theoretical results with the observations we perform a x 2 analysis, based 
on the expression 

2 



X 



E 



Mo.i Mo, 



(67) 



The quantities / xg ■ are the measured distance moduli for each of the supernovae of the 182 Gold06 



SNe la dataset 



771 ] . The jA i are the corresponding theoretical values and the a 2 i represent the 



measurement errors (cf. 



77 



78|). The probability density function (PDF) for our case with the 



free parameters Hq and qo is defined by 

P{H ,q ) = Ae- x2(Ho ' qo)/2 



(68) 



where A is a normalization constant. Marginalization over one of the two free parameters will 
lead to corresponding one-dimensional representations of the probability density. The details of 
the analysis are given in ref. [^sj]. The left panel of fig. [1] shows the two-dimensional probability 
distribution function for qo and h. Here h is defined by Hq = lOO/i km/Mpc/s. The best-fit value 
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FIG. 1: Left panel: two-dimensional probability distribution function for go and h. The brighter the color, the higher 
the probability. The best fit corresponds to go = —0.32 and h = 0.64. Center panel: one-dimensional probability 
distribution function for go. The probability is maximal for go = —0.39. Right panel: one-dimensional probability 
distribution function for h. The probability is maximal for h — 0.64. 

(minimum for x 2 ) is X 2 = 3.305, implying qo = —0.319 and h = 0.640. The corresponding values 
for the ACDM model are x 2 = 3.297 with f2 m o = 0.31, equivalent to qo = —0.535, and h = 0.645. 
For our power-law model, the present absolute value of the deceleration parameter is smaller than 
in the ACDM model. The one-dimensional probability distributions for q and h of our model are 
shown in the center and right panels, respectively, of figure[TJ The maximum values are qo = —0.392 
and h = 0.639. 

We conclude that even an approximate description of the dynamics on the basis of a power- 
law solution provides us with results that, at least on the background level, are consistent with 
observations. 



VII. DISCUSSION 



We have obtained exact power-law solutions for the dynamics of scalar-tensor theories for a 
constant ratio of the energy densities of the matter and the scalar-field components in the Einstein 
frame. Under this condition, the parameter that describes the interaction between matter and 
scalar field is no longer a free parameter. The Einstein-frame solutions are then transformed into 
the Jordan frame. The Hubble rates in both frames do not necessarily have the same sign. E.g., an 
expanding solution in the Jordan frame can correspond to a contracting solution in the Einstein 
frame. A corresponding property holds also for the deceleration parameter which may be negative in 
the Jordan frame while the corresponding quantity of the Einstein frame is positive. For pressureless 
matter the Jordan-frame EoS parameter w c s is related to its Einstein-frame counterpart w e g by 



18 



We fi = i_3jf ff • This solution defines the potentially possible effective EoS parameters for the 
cosmological dynamics as a whole. A de Sitter-type solution U7 e ff = —1 is obtained for u) e g = +1, 
which corresponds to a vanishing potential term in the Einstein frame. Other relevant cases, among 
them a matter-dominated phase, are recovered. The existence of exact power-law solutions with 
w e g < — ^ in a universe filled with pressureless matter can be seen as the simplest demonstration for 
the possibility to describe an accelerated expansion within scalar-tensor theories without a dark- 
energy component. We established an effective GR description of the Jordan- frame dynamics and 
identified the geometrical equivalent of dark energy. This component has a time dependent effective 
EoS, the time dependence being governed by w e s as well. The present value of the EoS parameter 
for such a dark-energy simulating component is consistent with the cosmological constant of the 
ACDM model. But the preferred present absolute value of the deceleration parameter is smaller 
than its ACDM counterpart. Our analysis is preliminary in the sense that it is restricted to the 
homogeneous and isotropic background dynamics. A perturbation analysis and a comparison with 
large-scale structure data will be the subject of future investigation. 

Acknowledgement: We are indebted to Julio Fabris for discussions and for providing us with 
the statistical analysis of Section IVII Financial support by FAPES and CNPq (Brazil) is gratefully 
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